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Abstract 

The study of laminated magnetostrictive functionally 
graded material (FGM) plate without shear deformation 
under thermal sinusoidal vibration and transient response 
is calculated by using the generalized differential 
quadrature (GDQ) method. In the thermoelastic 
stress -strain relations that containing a power-law function 
of a two-material FGM plate, the linear temperature rise 
and the magnetostrictive coupling terms with velocity 
feedback control. Four edges of rectangular laminated 
Terfenol-D FGM plate with simply supported boundary 
conditions are considered. The suitable product value of 
coil constant and control gain can be used to reduce the 
amplitude of center displacement into a smaller value. 
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Introduction 

Typical functionally graded material (FGM) is 
usually made of different phases constituent 
materials, for example, the ceramic and metal used in 
the engine combustion chamber to withstand 
ultra-high-temperature and to reduce the stress 
singularities, respectively. Chi and Chung (2006) 
presented the mechanical analysis of FGM plates 
subjected to transverse load. There were several 
studies in the transverse displacements for the 
laminated plate including the shear deformation 
effect. Amabili and Farhadi (2009) made the research 
of the shear deformable versus classical theories for 
nonlinear vibrations of rectangular isotropic and 
laminated composite plates. The over -prediction of 
natural frequencies in the solution without shear 
deformation was found. Ray and Shivakumar (2009) 
analyzed the effect of shear deformation on the 



piezoelectric fiber- reinforced composite plate by 
using the finite element method (FEM). Nguyen et al. 
(2008) obtained the static numerical results for the 
FGM plate with the effect of shear deformation. 

Magnetostrictive material Terfenol-D has the 
magneto- electric coupling property under the action 
of magnetism and mechanism. Hong (2009) used the 
computational generalized differential quadrature 
(GDQ) method to study the transient responses of 
magnetostrictive plates under thermal vibration. 
Thermal stresses and center displacement with and 
without shear effect were calculated in the thin and 
thick plate, respectively. Ramirez et al. (2006) 
obtained the free vibration solution for 
magneto-electro-elastic laminates through the Ritz 
approach. Lee and Reddy (2005) analyzed the 
non-linear response of laminated plate of 
magnetostrictive material under thermo-mechanical 
loading by using the FEM. Lee et al. (2004) obtained 
the transient vibration values of displacement for the 
Terfenol-D plate including the effect of shear 
deformation by using the FEM. Hong (2007) used the 
GDQ method to make the thermal vibration study 
for the Terfenol-D magnetostrictive laminated plate 
with the first-order shear deformation. Hong (2012) 
used the GDQ method to make the Terfenol-D FGM 
plate analyses under rapid heating induced vibration 
with the shear deformation effect. It is interesting to 
study thermal vibration of the transverse center 
displacement and thermal stress in the Terfenol-D 
FGM plate without the shear deformation effect by 
using the GDQ method. 

Formulation 

FGM 

Most materials of FGM can be used in the 
environment of higher temperature and can be 
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expressed in series form as follows by Chi and 
Chung (2006). 

i=\ 

(i) 

where P fgm is the material properties of FGM, n m 
is the number of materials mixed to form the FGM, 

V is the volume fractions, and > V- =1 for all 

i i 

constituent materials, P. is the individual 
constituent material properties, usually with the 
form as follows. 

p. ^p^p^r 1 + l + P x T + P 2 T 2 + P 3 T 3 ). 

(2) 

in which P ,P_ 1 ,P 1 ,P 2 and P 3 are the 

temperature coefficients, T is the temperature of 
environment. 

We denote the parameters for a two-material ( n m = 2) 
FGM plate as follows, a and b is the length in the 
x, y direction of the plate, n is the total thickness 
of magnetostrictive layer and FGM plate, h 3 is the 

thickness of magnetostrictive layer, h x and h 2 
are the thickness of FGM material 1 and FGM 
material 2, respectively, p x and p 2 are the in-plane 
distributed forces, q is the applied pressure load. 
The sum of volume fractions is in the form: 
V x +V 2 = 1 , the variation form of V 2 used in the 

power- law function is V 2 = (— ) Rn , where z 

h 

is the thickness coordinate, h is the thickness of 
FGM plate, R n is the power-law index. And the 
material properties for equation (1) can be assumed 
for the simple calculation and expressed as follows 
by Hong (2012). 

Efgm =(E 2 -E 1 )(l^f» + E l , 

(3a) 

v fgm =o 2 +n)/2, 

(3b) 

Pfgm =(P 2 +A)/ 2 ' 
(3c) 

a fgm =(«2 +«i)/2, 

(3d) 

shear stresses and shear strains in the laminate 



(3e) 

where E is the Young's modulus, v is the 
Poisson's ratio, p is the density, cc is the thermal 
expansion coefficients, k is the thermal conductivity, 
the subscript fgffl represents the FGM plate, the 
subscripts 1 and 2 represent the constituent material 
1 and 2, respectively. The property terms E x E 2 

V l V 2 A Pi a \ a 2 K \ K 2 are ex P resse d 

corresponding to term P. in equation (2). 
GDQ Method 

The GDQ method approximates the derivative of 
function, and the first-order and the second-order 
derivatives of function /*(jt, y) at coordinates 
(jc. , y j ) of grid point (/, j) can be discretized in 

series forms by Shu and Du (1997) and rewritten as 
follows: 



#1 

3c 



(4a) 



IV 1 



(4b) 



g I ~Y4 (2) f*, ■ 



1=1 



2 r* 



2 r * 



ckdy 



(4c) 

m=l 

(4d) 



/=1 



m=l 



(4e) 

where A. ( ^ and denote the weighting 

coefficients for the m th -order derivative of the 
function f*(x, y) with respect to x and y 
directions. 

Thermo Elastic Stress-Strain Relations with 
Magnetostrictive Effect 

We consider a rectangular laminated 
magnetostrictive FGM plate of the length a, b in the x, 
y direction, respectively, under uniformly 
distributed load and thermal effect. There are no 
without shear effect assumption. The plane stresses 
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in a laminated FGM plate with magnetostrictive 

effect for the k th layer are in the following 
equations by Lee and Reddy (2005): 
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(5) 

where a x and Ot y are the coefficients of thermal 
expansion, # is the coefficient of thermal shear, 
g. is the stiffness of magnetostrictive FGM plate, 
the simpler forms of Q. for FGM are given as 
follows. 



fin = fi 22 = 



J fgm 



l-v 



2 

fgm 

(6a) 



^ fgm E fgm 



1-v 



2 

fgm 

(6b) 



Q44 ~ Q55 ~ ~~ 



J fgm 



2d + v /s J 
(6c) 



and 



e 16 =e 2 6=fi4 5 =o. 

(6d) 

are the in-plane strains, 



AT = T (x,y,t) + —T l (x,y,t) is the temperature 
h 

difference between the FGM plate and curing area, 
1 is the coordinate in the thickness direction, h* is 
the plate total thickness, is the transformed 

magnetostrictive coupling modulus, f{ is the 
magnetic field intensity, expressed in the following 
equation by Lee et al. (2004), 

H z (x,yj) = kJ(x,yj) with velocity feedback 

control /(jc, y, t) = c(t) / in which k c is the c °il 

dt 

constant, Y(x, y,t) i s the coil current, c(t) is the 
control gain. 

Dynamic Equilibrium Differential Equations 

Without shear deformation effect, the time 
dependent of tangential displacement equations is 
assumed in the forms: u = u° (x, yj) 
v = v°(x,y,t) and the transverse displacement 
equation of the middle-plane is assumed in the form: 
W = w(x, y,t) in which t is time. The dynamic 
equilibrium differential equations in terms of 
displacements including the magnetostrictive loads 
are expressed in the following matrix forms by Hong 
(2009): 
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Where f 19 f 29 f3 are the expressions of thermal loads (TV, M) , mechanical loads p 2 , q) and 

magnetostrictive loads {N , M) , 



(7) 
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, dN x dN^ dN x dN^ 

fx = — - + + P\ + — - + ~ 

dx dy dx dy 



dx dy dx dy 



dx dy dy 

dM x 8M dM 
+ x - + 2 — + - 



dx dy dy 

(N x , M x ) = JJ + Q l2 a y + Q l6 a v )(T ,zj*T x )dz 

2 

i 

h* 

(N y ,M y ) = JJ (Q n a x + Q 22 a y + Q^a^ )(T ,zj^T x )dz 

2 

/ 

h*_ 

{N xy M^) = \^{Q l6 a x +Q 26 a y +Q^a xy ){T ,z^T l )dz 

2 

i 

h* 

(N x ,M x ) = jle 3l H z (l,z 2 )dz, 

2 

h* 

(N y ,M y ) = \^ 32 H z (l,z 2 )dz, 
(N X) ,,M xy ,) = fae 36 H z (l,z 2 )dz, 

2 

h* 

(A, , B tj ,Dy) = jj Q . (1, z , z 2 , (i , j = 1, 2, 6) , 

2 

/i* 

(/7, J ff) = jJ/? (l, Z )£fe. 

2 

in which p Q is the density of ply, p x and p 2 are 
the in-plane distributed forces, q is the applied 
pressure load. 

Dynamic Discretized Equations 

Without the shear deformation effect and under the 
following vibration of time sinusoidal displacement 
and temperature: 

u = [u°(x 9 y)]sin(a> mn t), 
(8a) 

v = [v°(x, v)]sin(^ mn 0/ 
(8b) 



w = w(x,y)sm(a) mn t), 
(8c) 

AT = [T (x,y) + ^T l (x,y)]sm(yt). 
h 

(8d) 

where the term co mn is the natural frequency of 
plate, y is the frequency of applied heat flux. 

We apply the weighting coefficients of discretized 
equations (4) in the two-dimensional GDQ method to 
discrete the differential equations (7). And we use the 
following non-dimensional parameters under 
sinusoidal temperature 

(7^ = 0, T x = T x sin(^- x/a)sin(^- y lb)). 

X -xl a, 
(9a) 

Y = y/b, 
(9b) 

U =u° I a, 
(9c) 

V = v°/b, 
_(9d) 

w = io/z*w/(^ x r 1 a 2 ). 

(9e) 

Thus, we can obtain the dynamic discretized matrix 
equations as in the paper by Hong (2009). 

Some Numerical Results and Discussions 

We like to consider the FGM plate consisting of two 
materials, of which the FGM material 1 is SUS304 
(Stainless Steel), the FGM material 2 is Si 3 N 4 

(Silicon Nitride). The temperature- dependent 
coefficients P Q9 P_ 19 P 19 P 2 and P 3 used to calculate 

material property terms, and E x E 2 V x V 2 p x 
p 2 GL X OC 2 and K x K 2 of these two typical 
constituent materials are listed in Table 1 by Shariyat 
(2008). The upper surface magnet ostrictive layer of 

the three-layer (0°™ / 0° /0 ) laminates FGM plate 
under four sides simply supported is considered, the 
superscript of m denotes magnetostrictive layer. 
The elastic modules, material conductivity and 
specific heat of the Terfenol-D magnetostrictive 
material are used the same value as in the paper by 
Hong (2007). We use the grid points for the GDQ 
computation as in the following coordinates: 
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x, = 0.5[1 - cos(-^— !- 7t)]a, i = l,2,-,N, 
N-l 

(10a) 

y ,. = 0.5[1 -cos{-^-7t)]b, = 1,2,...,M . 
; M-l 

(10b) 

The dynamic convergence results are obtained for 
center displacement amplitude w {a 1 2, b 1 2) 
without shear effects in the thermal vibration of 
sinusoidal temperature only (T = 0, T x =100°^, 
p x = p 2 = q = ) at time t = 6s, mode shape 
m — n — \ , with control gain value k c c(t) = , 
aspect ratio alb — 0.5, 1 and 2, side-to-thickness 
ratio alll = 100, 50, 20, 10 and 5, /** =1.2 mm, 
/z 3 = 1 mm, h ] = h 2 = 0.1 mm, R n - 1 , 

7 = 300°^. Table 2 shows the w(a/2,b/2) 
(unit mm) in the grid point NxM = 17x17 

21x21 25x25 29x29 and 33x33 of GDQ 
method for the laminated Terfenol-D FGM plate at 
time t = 6 s. We find the NxM = 33x33 grid 
point have the good w (a/2,b/2) convergence 
result and use further in the GDQ analyses of time 

responses for displacement and stress for a! h = 5, 
20 and 100. 

We obtain the lowest frequency y of applied heat 
flux and vibration frequency 0) n of Terfenol-D 
FGM plate ( m = n = 1 ), at time t = 0.001s, Is, 
2s,. ..and 9s, NxM = 33x33, h* = 1.2 mm, 
h 3 = 1 mm, h ] = h 2 = 0.1 mm, alb — \ , q = 0, 
R n = \, T = 653°K, f x =100° K, as shown in 
Table 3. For the magnetic coil constant k c effect on 
the w {a 1 2, b 1 2) (unit mm) of vibration under the 
constant gain value c(t)= 1, time t= 6 s, a I h =5 
without shear, we obtain the sketch of w(a 1 2, b 1 2) 
vs. k c as shown in Fig. 1. The suitable product 
value of k c and c(t) can be used to reduce the 
amplitude of w (a 1 2,b 1 2) into a smaller value 
near 0.00, there is an amplitude peak value at k c c{t) 

= 7.16E08, we find k c c(t)= 1.56E09 for al\l =5 
thick plate can be chosen as the best suitable product 
values. 

Firstly, thermal sinusoidal vibration is investigated 
with time step equal to 0.1s, the suitable chosen 
product k c c{t) of coil constant and controlled gain 



values versus time t for R n =1 , thick plate 
a l}{ =5 and thin plate a/h*= 20,100 at 
T = 653 K, as shown in Table 4. 

Fig. 2 shows the w (a/2,b/2) (unit mm) versus 
time t of GDQ method for the laminated Terfenol-D 

FGM plate a III =5 , 20 and 100, respectively 
without shear effects. At time t= 0.001s, there is a 
great amplitude value of displacement with 
uncontrolled value ( k c c(t) = 0), w (a/2,b/2) = 

-0.156429 mm for thick plate a/h*=5 , 
w(a/2,b/2) = -4.31174 mm for thin plate 

a III =20, w(a/2,b/2)= -022975 mm for thin 

plate a III =100 . We find the amplitudes of 
W (a/2,b/2) with controlled k c c(t) values are 

smaller than the amplitudes of w (a / 2,b / 2) with 

uncontrolled value (k c c(t)= 0), generally by using 

the GDQ method. We can use the suitable product 
values of controlled k c c(t) to reduce the amplitude 

of w (a 1 2, b 1 2) into a smaller value near 0.00. 

Fig. 3 shows the time response of the dominated 
dimensional stress <j x (unit GPa) at center position 

of upper surface Z = 0.5/z with respective to time 
for the laminated Terfenol-D FGM plate a III =5, 
20 and 100, respectively without shear effects 
a =0 and <j — . We find the maximum 

xy xy 

response values of <j x = -9.13E-04GPa for a I h* = 5 
and 20, a x = -9.12E-04GPa for a//z*=100 with 
controlled k c c(t) case are almost equal to the 
response values of <j x = -9.12E-04GPa with 
uncontrolled case ( k c c(t) = 0), generally by using the 
GDQ method. 

Fig. 4 shows the compared non-dimensional 
W(X,b/2,6) versus X of a III =5, 20 and 100 
without / with shear effects, for Terfenol-D FGM 
plate, at time t= 6s, m = n = 1 , NxM= 33x33 
for the case of without shear and 17x17 for the 
case of with shear, II =1.2 mm, h 3 =1 mm, h x = 
h 2 = 0.1 mm, alb = 1 , q = , R n =1 , 
T = 653° K, T x = 100° K, k c c(t) = 0. In the shear 

effect case, we use the YNS first-order shear 
deformation theory for the time dependent of 
displacement field, the value for shear correction 
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coefficients k a = — 516 used in the dynamic 
equilibrium differential equations by Hong (2012). 
For the thick plate a I h = 5 , the values of 
displacement W(X,b/2, 6) versus X without 
shear case are smaller than that with shear case. The 
maximum value W(X,b/2,6) = 0.0821655 occurs at 

X = 0.5 of the with shear case, and W(X,b/2,6) = 
-0.0325304 occurs at X = 0.691342 of the without 

shear case. For the thin plate a I h = 20, the values 
of displacement W(X ,b/2,6) versus X without 
shear case are greater than that with shear case. The 
maximum value W(X,b/2,6)= -0.00133522 occurs at 

X = 0.777785 of the without shear case, meanwhile 
W(X,b/2,6)= 0.00033884 occurs at 1= 0.5 of the 

with shear case. For the thin plate a I h =100, the 
values of displacement W(X,b/2, 6) versus X 
without shear case are greater than that with shear 
case. The maximum value W(X,b/2,6) 
0.0000138781 occurs at 1 = 0.0380602 of the without 
shear case, while W(X,b/2,6) = 0.000000541036 
occurs at X = 0.5 of the with shear case. The values 
of displacement W(X,b/2, 6) are decreasing with 

a I h values increasing for square plate alb = 1 
with and without shear deformation effects. The 

deflection of center position versus alh at t= 6 
sec in Fig. 4d shows the deflection of thick plate 

alh =5, and the case of deflection value with 
shear effect is much greater than the case of without 
shear effect. The deflections with and without shear 
effects are almost in the same values for the thin plate 

alh" =100. 

Secondly, transient response is investigated with time 
step equal to 0.001s and use the fixed frequency 
/ = 523.599/s of applied heat flux. The suitable 
chosen product k c c{t) of coil constant and 
controlled gain values versus time t for R =1 , 

thick plate alh* =5 and thin plate alh* =100 
at T = 653° K , as shown in Table 5. We find the 
more thin plate is easier to control the displacement 
with less values of k c c{t) . 

Fig. 5 shows the transient value w (a 1 2,b 1 2) (unit 
mm) versus time t of GDQ method for the 

laminated Terfenol-D FGM plate a I h* =5 and 100, 
respectively without shear effects. At time t= 0.001s, 
there is a great amplitude value of displacement with 



uncontrolled value ( k c c(t) = 0), from w (a 1 2, b 1 2) 
= -0.156429 mm converges to small value for thick 
platea//** =5, from w (a/2,b/2) = -0.22975 mm 

converges to small value for thin plated/ h = 100. 
We find the amplitudes of W (a/2,b/2) with 

controlled k c c(t) values are smaller than the 
amplitudes of w(a/2,b/2) with uncontrolled 
value ( k c c(t) = 0), generally by using the GDQ 
method. We can use the suitable product values of 
controlled k c c(t) to reduce the amplitude of 
W (a 1 2, b 1 2) into a smaller value near 0.00. 

Fig. 6 shows the transient value of the dominated 
dimensional stress <y x (unit GPa) at center position 

of upper surface Z = 0.5/z versus time for the 
laminated Terfenol-D FGM plate a I h* =5 and 100, 
respectively without shear effects = and 

<7 xy = . We find the transient values <j x oscillate 
between -9.12E-04GPa and 9.12E-04Gpa for 
alh —5 and 100. The transient values <y with 
controlled k c c(t) case are almost equal to the 
values with uncontrolled case ( k c c(t) = 0), generally 
by using the GDQ method. 

Conclusions 

The GDQ calculation provides a method to compute 
the controlled displacement and stress in the 

(0° m /0°/0°) ply Terfenol-D FGM plate subjected 
to thermal vibration and transient response of 
sinusoidal temperature without shear deformation 
effect. The computation provides the following 
results, (a) The suitable controlled product values of 
coil constant and control gain k c c(t) can be used to 

reduce the amplitude of center displacement 
W (a 1 2,b 1 2) into a smaller value near zero, (b) 
The amplitudes of stresses <y x of plates are almost 

in the same values under the two cases: with and 
without k c c(t) values, (c) The values of 

non-dimensional displacement W(X,b/2, 6) are 

decreasing with alh values increasing for square 
plate a/b = l under the two cases: with and 
without shear deformation effects, (d) The deflection 

values of center position versus a/h* at t= 6 sec 
are investigated, which found that the value with 
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shear effect is much greater than the case of without 

shear effect for the thick plate a I h = 5 , as well the 
deflections with and without shear effects are almost 

in the same values for the thin plate a I h = 100 . 
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FIG. 2a w(a/2,b/2) vs. t for dill = 5 



FIG. 3a <T vs. t for a I ' h 
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fig. 2b w(a/2,b/2) vs. t for alh* = 20 
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FIG. 3b <y vs. t tor alh = 20 



0.05 r 



0.00 
-0.05 
-0.10 
-0.15 
-0.20 
-0.25 



* uncontrolled (* c c(0 = 0) 
— controlled (with MO) 



0.00E+00 
-1.00E-04 
-2.00E-04 
-3.00E-04 
-4.00E-04 
-5.00E-04 
-6.00E-04 
-7.00E-04 
-8.00E-04 
-9.00E-04 
-1.00E-03 



uncontrolled (k c c(t) = 0) 
controlled (withMO ) 



1 2 3 4 5 6 7 



0123456789 



FIG. 2c W (a/2,b/2) vs. t for alh =100 

fig. 2 w(al2,bl2) vs. t foxalh* =5,20andioo 

WITHOUT SHEAR 



FIG. 3c <J x vs. t for alh =100 

FIG. 3 CT X vs. t foxdlK = 5, 20 AND 100 WITHOUT 
SHEAR 
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W(X,b/2, 6) 




0.0 0.2 0.4 0.6 0.8 



1.0 



FIG. 4a W(X,b/2,6) vs. X for a I ' h* = 5 
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_aJh 
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FIG. 4d CENTERDEFLECTION vs. alh AT f=6sec 

FIG. 4 COMPARED W(X,b/2, 6) vs. X iorCllh* =5, 
20 AND 100 
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FIG. 4b W(X,b/2,6) vs. X for =20 
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FIG. 4c W(X,b/2,6) vs. X for alh* =100 



FIG. 5a TRANSIENT VALUE OF w{d I 2,b I 2) vs. t FOR 

alh* =5 
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FIG. 5b TRANSIENT VALUE OF w(a/2,b/2) vs. t for 

alh* =100 

FIG. 5 TRANSIENT VALUE OF w(a/2,b/2) vs. / for 
a / h* =5 AND 100 WITHOUT SHEAR 
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■ k c c{i) = 
- with k c c(t) 




111 1,111 11 
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FIG. 6a TRANSIENT VALUE OF <j vs. / FOR a/h* =5 




t,cffl = 

— - — with k c c(t) 




III 
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FIG. 6b TRANSIENT VALUE OF <J x vs. t FOR 

a/h* =100 

FIG. 6 TRANSIENT VALUE OF <J x vs. t FOR a/h* 
AND 100 WITHOUT SHEAR 



TABLE 1 TEMPERATURE-DEPENDENT COEFFICIENT OF CONSTITUENT FGM MATERIALS 



Materials 




P 


Pi 


Pi 


p 2 


p 3 


SUS304 




201.04E09 





3.079E-04 


-6.534E-07 









0.3262 





-2.002E-04 


3.797E-07 







p^Kg/m 3 ) 


8166 
















a x (°K~ l ) 


12.33E-06 





8.086E-04 










k x {W lm°K) 


15.379 
















C vl (JIKg°K) 


496.56 





-1.151E-03 


1.636E-06 


-5.863E-10 


Si 3 N 4 


E 2 (P fl ) 


348.43E09 





-3.70E-04 


2.16E-07 


-8.946E-11 




^2 


024 
















p 2 (Kglm 3 ) 


2370 
















a 2 CK' 1 ) 


5.8723E-06 





9.095E-04 










k 2 (W lm°K) 


13.723 
















C v2 (J/Kg°K) 


555.11 





1.016E-03 


2.92E-07 


-1.67E-10 
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TABLE 2 DYNAMIC CONVERGENCE OF TERFENOL-D FGM PLATE WITHOUT SHEAR 


alh* 


GDQ method 




W {a 1 Z, u 1 L) (unit mm) 




NxM 




alb = 0.5 


alb = 1.0 


alb = 2.0 


1 nn 
1UU 


17x17 




-U.oo41oOe.-U4 


n nn OT7 c\a 
-U . y\J 04 1 ZE -U4 


n 1 1 K£H o~c no 
-U.l/OO/oE-Uo 




21x21 




n caisson r\A 

-U.5o3ooUh-U4 


n 'vn/^oi nc n/i 

-U./9oolUh-U4 


n lA/i/'TTC no 

-U.lo4oz/E-U3 




25x25 




n d oir\ r TC n/i 

-U.o1o1U/e-U4 


n 'vo d r 7ni7 n/i 

-U./Jol/UL-U4 


n 1 non 1 17 no 
-U.loUz9lE-UJ 




29x29 




-U.49459ZE-U4 


n 'vnn'V/i ci7 n/i 
-U./U9/45E-U4 


n 1 c eve nc no 
-U.Ioo/oUe-UJ 




33x33 




n zioooonc r\A 

-U.4od3zUh-U4 


n r 7nio r 7cc n/i 

-U./U 1J/5L-U4 


n 1 c ooo nc no 

-U.15oZoUe-U3 


en 


17x17 




n OnOTI OT7 HQ 

-U.zUz/Ioe-Uo 


n OQ001 of no 
-U.zooolZE-Uo 


n cq ceo ^c no 
-U.OOOOOOE-UJ 




21x21 




n ini ct/T no 


n o'vi s aitj no 
-U.Z/ 1641L-UJ 


n C7mi nc no 
-U.o/9z1Ue-Uv3 




25x25 




-U.lo/OooE-Uo 


-U . ZO yD\J Ze -Uo 


n C77CQaT7_no. 




29x29 




-U.loOOOOE-UO 


n o^qi /i nc no 
-U.Zo^MUe-Uo 


n c77 , 2Qor no 

-U.3/ /o>>ZE-Uv3 




33x33 




n lo/m r\jj no 
-U.lool9Uh-UJ 


n o/^nci 1 T7 no 
-U.zo9ollE-Uv3 


n cu no 
-U.o//4zoE-Uv3 


20 


17x17 




-0. 99911 6E -03 


r\ 1 o 1 no /r no 

-0.131836E-02 


r» rnn or no 

-0.205093E-02 




21x21 




-0.987938E-03 


-0.131330E-02 


-0.20 505 OE-02 




25x25 




n ooiC/i/^T7 no 

-U.9o1o4oe-U3 


n 10100 0T7 no 
-U.L3 UoJe-Uz 


n oncn/i /ic no 
-U.zUoU44e-Uz 




29x29 




r\ nnn/io/^c no 

-U.99U430E-U3 


n ioin/iof no 
-U.L3 1U4ZE-UZ 


n on enn /i c no 

-U.ZU5U94E-U2 




33x33 




n nn c\~!C cwn no 
-U.99U/oUh-UJ 


n 1 o nno q~c no 
-U.UUUOOE-Uz 


n on cm nc no 
-U.ZUoUIUe-Uz 


1 n 


17x17 




-U.Z41Z0OE-UZ 


_n OQ^o^T7_no 
-U.ZoooZoe-UZ 


n OQo.Qo.ni7 no 
-U.Z^oooUE-UZ 




21x21 




n O/i 1 oq /if no 
-U.z4L3o4h-Uz 


n OQoenne no 
-U.zoJoUyE-Uz 


n onoococ no 
-U.zyjoooE-Uz 




25x25 




_n 9A 1 1 r 1 no 

-U .Z4 1 ID 1E-UZ 


-U . Zo OOO OE -UZ 


-U . Ay 0040E-UZ 




29x29 




-0.241165E-02 


-0.283456E-02 


-0.293842E-02 




33x33 




-0.240542E-02 


-0. 28341 OE-02 


-0.294010E-02 


5 


17x17 




-0.340208E-02 


-0.343168E-02 


-0.232384E-02 




21x21 




-0.340130E-02 


-0.342991E-02 


-0.232400E-02 




25x25 




-0.340583E-02 


-0.343044E-02 


-0.232387E-02 




29x29 




-0.340369E-02 


-0.343035E-02 


-0.232335E-02 




33x33 




-0.339761E-02 


-0.343055E-02 


-0.232343E-02 




TABLE 3 y 


AND 0) n OF TERFENOL-D FGM PLATE WITHOUT SHEAR 




a/h* 


f =0.001 s 


t = 


Is 


t=2s 


r 


co n 


r 


r 




1 nn 
20 

5 


523.599 
523.599 
523.599 


0.569829E-02 
0.567729E-02 
0.196036E-01 


1.57080 
1.57080 
1.57080 


0.141884E-01 0.785398 
0.151387E-01 0.785398 
0.136018E-01 0.785400 


U.J070JUCDZ 

0.567730E-02 
0.5381 13E-02 


alh* 


t = 


3s 


t = 


4s 


t=5s 


r 


co n 


r 


r 


co n 


100 
20 
5 


0.523599 
0.523599 
0.523599 


0.569830E-02 
0.567730E-02 
0.538112E-02 


0.392699 
0.392699 
0.392701 


0.569830E-02 0.314159 
0.567730E-02 0.314159 
0.538113E-02 0.314161 


0.569829E-02 
0.567729E-02 
0.5381 13E-02 


alh* 


t = 


6s 


t = 


7s 


?=8s 


r 


co n 


r 


r 


^11 


1 nn 
20 

5 


0261799 
0.261799 
0.261800 


0.569830E-02 
0.567730E-02 
0.5381 13E-02 


0.224400 
0.224400 
0.224401 


0.569830E-02 0.196349 
0.567730E-02 0.196349 
0.538112E-02 0.196351 


n R^QS^nF-n? 

U.J070JUL-UZ 

0.567730E-02 
0.5381 13E-02 


alh* 


t = 

r 


9s 








100 
20 

5 


0.174532 
0.174532 
0.174534 


0.569830E-02 
0.567730E-02 
0.5381 13E-02 
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TABLE 4 k c C(t) VERSUS £ FOR R n =l WITHOUT SHEAR 



al\i 






£ e c(0 








£ = 0.001 -0.1 


t=02~0A 


? = 0.5 -0.9 


t = 1.0 t = 


1.1-1.4 £ = 15 


5 


-10 9 


-3xl0 6 


-10 7 


-10 2 


-10 9 -10 7 




£ = 1.6- 3.3 


? = 3.4- 3.9 


£ = 4.0- 4.6 


£ = 4.7- 9.0 






-10 9 


-2xl0 9 


10 9 


1.56xl0 9 






£ = 0.001 


t = 0.1 t = 0.2 


£ = 0.3-0.9 


£ = 1.0 t=ll 


£=1.2-1.3 




-10 9 


10 6 -10 9 


-10 5 


-10 2 ~10 9 


-10 7 




£ = 1.4- 3.2 


? = 3.3~ 3.5 


? = 3.6 


£ = 3.7- 5.2 


t = 5.3- 5.4 


20 - 


10 9 


10 5 


10 2 


-10 9 


10 5 




£ = 5.5- 7.5 


t = 7.6~ 9.0 










-10 7 


10 9 










£ = 0.001 


t = 0.1-1.5 


1.6 -1.9 


£ = 2.0 £ = 2.1- 2.3 


/" = 2.4 


100 


-10 9 


10 2 


-10 9 


-10 2 10 7 


-10 9 



£ = 2.5-9.0 
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